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Abstract. We study the effect of resonant scatterers on the local density of 
states in a rectangular graphene setup with metallic leads. We find that the 
density of states in a vicinity of the Dirac point acquires a strong position 
dependence due to both metallic proximity effect and impurity scattering. This 
effect may prevent uniform gating of weakly-doped samples. We also demonstrate 
that even a single-atom impurity may essentially alter electronic states at low- 
doping on distances of the order of the sample size from the impurity. 
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1. Introduction 

Since the discovery of graphene [H E] a great progress has been made in production of 
high quality samples with nearly ballistic behavior of Dirac electrons O HI [5] . Even 
though the experimental results agree well with the ballistic theory of transport [5] [7] , 
the role of metallic contacts in electronic properties of sub-micron graphene samples 
is debated [8]. In this paper we introduce a notion of metallic proximity effect in 
graphene and show that this effect can lead to a non-uniform gating of small graphene 
flakes. 

Despite the intense graphene studies, the scattering mechanisms, which limit 
electron mobility in the best available graphene flakes, are not clearly identified [S]. 
Still, the experimental evidence that strong impurities, such as adatoms and vacancies, 
dominate the scattering at low temperatures is currently building up [TO] EJ Q21 Q2] . 
We demonstrate below that even the point-like impurities may introduce large spacial 
variations of electronic density at distances of the order of the system size from the 
impurity provided the quasiparticlc energies are close to the Dirac degeneracy point. 

An impurity can be characterized by at least two length scales: its physical size a 
and the characteristic scattering length l s . The latter is determined by the scattering 
cross-section and, therefore, depends on the energy of Dirac quasiparticles that scatter 
on the impurity. A single weak impurity, i.e. the impurity with l s -C a, has practically 
no effect on transport and electronic density of states in ballistic sample. The effect 
becomes essential only for a substantial concentration of weak impurities [TH [T71 [TS] . 
In contrast, a single strong impurity, i.e. the one with £ s 3> a, can have a noticeable 
effect on transport [15] . 

The strong impurities are of a particular importance at zero doping, i.e. in the 
situation when the chemical potential in graphene is tuned precisely to the Dirac point. 
It is theoretically predicted that each strong impurity can form a resonant "mid-gap" 
state [14], which contributes a universal value to the conductance and shot noise [15] 
provided the wave-length of the Dirac quasiparticle, A, exceeds the distance between 
the metallic leads, L. The resonant conditions are equivalent to the limit £ s 3> L. 
Evidently, in this regime, the properties of graphene flakes are strongly influenced by 
both the impurities and the sample boundaries. 

The values of conductance and shot noise observed in experiments with sub- 
micron graphene flakes [4] [5] agree with the predictions of ballistic theory [6] [7] for 
A > L, which makes us to believe that this limit is experimentally achievable. Both 
the ballistic theory of ref. [7] and experiment [H [5] deal with rectangular graphene 
samples with large aspect ratio, W/L 3> 1, where W is the sample width and L is the 
distance between the metallic leads. 

For A ^ L the proximity to metal gives rise to the so-called pseudo-diffusive 
charge transport [19] in the model with abrupt boundary conditions at graphenc- 
metal boundaries. We shall demonstrate below that in this model the local density 
of states (LDOS) takes on the most simple form, po(x) = 2[irhvL sin(7rx / 'L)]^ 1 , which 
corresponds to the lowest curve in figure [1] The enhancement of the density near the 
graphene-metal interfaces illustrates the metallic proximity effect. The LDOS in the 
presence of a single impurity is shown in figure [2] with a contour plot for resonance 
and off-resonance conditions. Such a strong spacial dependence of the LDOS suggests 
that, in the presence of electron-electron interactions, the sample is not uniformly 
charged by the gate electrode. 
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2. Clean sample 

The LDOS is defined by the general formula, p(r;s) = — 7T _1 Im tr G R (r, r; e), where 
the retarded Green's function yields the matrix equation 

(e + iO - n(x) - H)G R (r, r') = S(r - r'), (1) 

with H = Hq + U(r). In this paper the Dirac Hamiltonian of ideal graphene, 
Hq = —ikvcrV , is taken in the valley symmetric representation, where v = 10 6 m/s 
is an effective "speed of light". The impurity potential, U(r), is considered to be a 
scalar in the Dirac equation, which assumes that the corresponding physical potential 
is smooth on the scale of the graphene lattice. The vector <x = (a xl o~ y ) is a two- 
component vector of Pauli matrices and the energy, s, determines the wave-length of 
free Dirac quasiparticle as \e\ — hv/X. 

Following earlier works we describe the metal-graphene interfaces at x = 0, L with 
a step function model 

, , _ / -hv/Xp, x < 0, x > L . . 

^ x >-\ 0, 0<x<L ' (Z) 

where Xf is regarded as the Fermi wave length in the metallic lead. We naturally 
assume that Xf <C min{L,A}. In order to avoid a detailed discussion of the sample 
edges we impose periodic boundary conditions in y direction. 

We start with a solution of equation |T|) in the absence of impurity potential. 
This solution is most easily obtained in the channel representation, which is defined 
by the Fourier transform, 

G R (r, r') = ^ Yl e tq(y ~ y,) G R (x, x'), (3) 

where the summation runs over the discrete values of transversal momentum, q n = 
2im/W . Selecting the solutions which decay in the leads we obtain the effective 
boundary conditions on graphene-metal interfaces in the limit Xp <C min{L, A} |15j . 

(1,1^(0,0=0, (l,-l)Gf(L,x')=0, 

G«(:r,0)(l,-1) T = 0, G R ( Xl L)(l,if = 0. (4) 

Using these boundary conditions the Green's function inside the clean sample is readily 
found as 



R , , , _ -i(F 1 (u 1 )+F 2 (u 2 )) 
^ , q (x,x,e)- 2hvcos9cos $ ■ 

where u\ = $ — k\x — x'\, u 2 — k(L — x — x') and 



(5) 



k = y/{e-i0) 2 - (hvq) 2 , e lB = hv(k + iq)/e, (6) 

e « =e <W(l_ e *«)/(l + e <fl ). (7) 

In equation (|5|) the following definitions were used, 

F( u )=( -isinu rj cos(u + r] 9) \ , g . 

1 ^ ' I r\ cos(u — r\ 6) — isinu J' ^ ' 

F 2 (u)= ( -- in (- + ) -«»« V ( 9 ) 

v ' V cos it ism{u — 9) J K ' 

where r\ = sign (x — x'). The quantity, k, introduced in equation ([5]) has a meaning of 
x component of the quasiparticle momentum, which can take on complex values. 
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Thus, the LDOS of a clean system is given by the expression 




(10) 



where the factor 4 takes into account the spin and valley degeneracies. The LDOS 
of the clean setup has no dependence on y due to the translational symmetry of the 
system in this direction. The x-dependence of the LDOS is plotted in figure [T] for 
W/L = 10 and several values of energy. 

The expression (|10p can be approximated in the doped graphene, A <C L, by 
averaging over the Friedel oscillations on the scale A. Away from the boundaries, 
d = min{a;, L — x} A, the evanescent modes, which are characterized by complex 
momenta k, can be disregarded and the averaging over the Friedel oscillations can 
be replaced by the averaging over the phase <j> = kL. For W ^> L, the averaged 
LDOS value is given by po = 2\e\/(irhv) 2 , which is shown in figure [1] with dotted 
lines. Near the sample edge, d = minja:, L - i} < A, the density is dominated by the 
metallic proximity effect, so that it diverges at the metal-graphene interfaces in the 
limit Xp — ► 0. For d <C A = hv/\e\ and W ^> L, the LDOS does not depend on energy 
and can be approximated as po — 2[ir 2 hvd}~ 1 . 

In contrast, for X ^S> L, the LDOS is fully determined by the evanescent modes. 
In particular, for e = (A — > oo), the result can be obtained by taking k = —iq, 
9 = +ioo, and <I> = —iqL in equation (|10p . which leads to the expression 



In the limit W 3> L, the sum over the transversal momenta, q, can be replaced by the 
integration, with the result 



where the small term A^/4 in the denominator is effective only at the edges, x = 
and x = L, in which case the sum (jlip equals the number of channels in the lead. 

The result (fl"2"|) shows that the LDOS at e = is dominated, in the limit W ^> L, 
by the states originating in the leads. These boundary states, which are also referred 
to as the evanescent modes, are responsible for the so-called pseudo-diffusive transport 
regime in this setup |19j . which is characterized by universal full counting statistics [15j. 
We see, therefore, that the strong metallic proximity effect is the natural property of 
the abrupt interface model ([2|), which is compatible with the transport measurements 
of ref. 0. 

3. Effect of impurity 

In this Section we consider the effect of a scalar point-like impurity on the LDOS at 
e = for a rectangular graphene sample with W ^> L. The calculation of the Green's 
function is analogous to the one presented in ref. |15j for the charge transport and can 
be carried out straightforwardly. 

We shall start with the real space representation of the free Green's function ([5|) 
at e = 0, that reduces, in the limit W ^S> L, to [15] 




(11) 



Po(x) 



2 



(12) 



irhv(Lsm(nx/L) + Af/4) ' 




(13) 
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C±(r,r') = [Lsw.{n{x±x' + i{y-y'))/2L)+\ F /A]~ 1 . (14) 

One can see from the equations (|13I14|) that the free Green's function obeys the chiral 
symmetry 

[G*(r',r)f = -a z G«(r,r<)a z = G£(r,r'), (15) 

which is, therefore, preserved in the abrupt interface model described by equation 1(2]). 

One can also see from equation (fT5]) that the Green's function at the coinciding 
arguments can be decomposed into the sum of the free Green's function in an infinite 
system, g, and the proximity-induced part, G rcg , such that 

km GX(r,r')=g(r,r') + G Ies (r), (16) 

r' — >r 

where g is off-diagonal and singular, 

9(r,r ') = -— — ^ 17 

2nhv \r — r'Y 

while G^ g (r) is proportional to the product of the LDOS (fl"2]) and the unit matrix, 
Gf og (r) = (-m/4)p (x). (18) 

In order to simplify the calculation we choose the impurity potential in the form 
of a disk with sharp boundaries, U(r) = u 9(a — \r — ro\), where the theta function, 
0(x), equals 1 for x > and zero otherwise. We also take the limit of point-like 
impurity, a/L < 1, In this limit the full T-matrix of the impurity can be written, in 
the symbolic notations, as 

T(r ) = U - 1 7 ^T T =T — T . 1 , (19) 

1 - G$U 1 + («7r/4) j o (xo)T 



where we used the decomposition (|16[) and introduced the reduced T-matrix, T = 
U[\ — gU]^ 1 , which is represented, in the detailed notations, by the infinite sum 

oo r s 

T = J2 H^rpUir^gir^^Uirz) . . . g(r s -i,r s )U(r s ). (20) 

8=1"' p=l 

The reduced T-matrix given by equation (|20p does not contain any information 
on the metallic leads and sample boundaries. Therefore, it can be obtained from the 
consideration of scattering of free Dirac quasiparticles on a single disk impurity. This 
problem has been addressed by several authors [20] [211 US • The result in the s-wave 
approximation and in the limit e —> takes the form 

T =2*1*1., ^8 = 4^, ( 21 ) 

Jo(ua/hv) 

where J n is the Bessel function and the T-matrix, T, is proportional to the unit matrix 
in spin, valley, and sublattice space. 

The equation (pH) completes the calculation of the Green's function in the 
presence of a single impurity, 

G R (r,r') = G«ir,r')+G*ir,r )T(r )GX(r ,r), (22) 

where r = {x ,y ) specifies the impurity location. We stress that the characteristic 
impurity size, t s , determined by the reduced T-matrix, strongly depends on the 
impurity potential, u, and diverges at the resonance. 
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From equation (|22|) we find the LDOS at the point r, which can be written as 

n n , , JC_(r,r )| 2 - \C+{r,r )\ 2 

p(r)=p (x)+po(xo) 5-. (23) 

(2/tt4) 2 + [7rHvp (xo)/2} 2 

This result is plotted in the figure[5]for an impurity placed in the middle of the sample, 
x Q = L/2. 

For the resonant conditions, J {ua/hv) = 0, (i.e. for £ s — » oo) the disk impurity 
supports a quasi-bound state at e = [551 EH] with only one pseudo-spin component 
being localized. In this case, the LDOS is strongly enhanced at distances of the order 
of L from the impurity. Away from the resonant condition, the effect of a single 
impurity on the LDOS is predictably small. It was argued on the basis of the DFT 
analysis [14j that some single-atom impurities such as hydrogen naturally create the 
mid-gap quasi-bound states. The resonance condition for such scatterers are always 
fulfilled near the Dirac point. Therefore, even few impurities in the sample will have a 
great impact on the LDOS. We expect that this effect can be observed with the STM 
measurements. 

Finally we note that the denominator in the second term of equation (|2"3"]) coincides 
with the denominator in the impurity contribution to the conductance in the same 
setup [T5]. Near the resonance, the second term in equation (f2"3"]) has a Lorentzian 
dependence on the position of the impurity level with respect to the Dirac point. 
This characteristic dependence reflects the resonant nature of charge transport that 
originates in the impurity assisted tunneling in undoped graphene 23J. 

4. Discussion 

We found from equation (p~0|) that the LDOS in ballistic rectangular sample with 
W ^> L can be approximated as 

2 f hv' 

where d = min{a;, L — x} is the distance to the closest graphene-metal interface. This 
result shows that the metal proximity effect is fully dominating the LDOS for energies 
close to the Dirac point, |e| < Hv/L. It also demonstrates that the LDOS is never 
vanishing in graphene with metal boundaries. 

From equation ([24]) we readily estimate the integrated density of states as 



p a (x; e) = ,_ fc , 2 max \ \e\, — \ , (24) 



2W 

»° £ =-a»r + 2 1h(£/Af )), (25) 

where £ = min{L, 7rftw/|e|}. Thus, for |e| < Hv/L, the integrated density increases 
logarithmically with the system size, L, due to the electronic states at the metal- 
graphene interfaces. Peculiar transport properties of ballistic graphene [HI [15] are 
originating in the metallic proximity effect. 

Similarly, the presence of the resonant impurity has a great effect on the LDOS 
near the Dirac point since the corresponding electronic states localize not only on the 
metal interfaces but also on the impurity. This effect leads to a highly inhomogeneous 
electronic density profile for |g| < Hv/L, which is illustrated in figureO 

This result implies that electron-electron interactions may prevent the uniform 
gating of the sample in the vicinity of the Dirac point. Assuming that the interaction 
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potential is given by V(r) we can estimate, up to a constant value, an additional 
potential profile in the single-electron problem ([T]) as 



where /(e) is the Fermi distribution function. Since only the spatially varying part 
of the potential <f>(r) is relevant, we can let E c = Hv/L. Finally the LDOS and 
the potential profile, <fi(r), have to be found self-consistently from equation |T]) with 
/J>(x) — > n(x) + <j)(r) and equation (|26[) . The detailed analysis of the non-uniform 
gating is, however, beyond the scope of the present paper. 
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0.0 0.2 0.4 0.6 0.8 1.0 
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Figure 1. The local density of states of Dirac quasiparticles from the equation 
lllOI I in ballistic graphcne sample with metal boundaries and W/L = 10. Different 
curves correspond to different energy values eL /hv = 0, 2, 4, 8. In the middle of the 
sample (i.e. for J«i<L-A, where A = 2irhv/\e ) the LDOS is approximately 
given by po = MA 71 "^) 2 (dotted lines). Near the metal boundaries (with \p = 0) 
the LDOS behaves as po = [it 2 hvd\~ 1 , where d = min{x, L — x}. 
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Figure 2. The density plot for the LDOS at e = calculated from the equation 
J23} or out-of-resonance (left) and resonance (right) conditions. The physical size 
of impurity is chosen as a = 0.05L, while the effect on the LDOS is determined 
by the effective length £ s . At resonance conditions, i s 3> L, the LDOS is strongly 
enhanced in the large area, which is of the order of the system size L. 



